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Precision and accuracy of the orbital parameters derived from 
2D & 1 D space observations of visual or astrometric binaries 
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Abstract. Recent investigations (e.g Han et al. 2001 ) have shown that fitting the Hipparcos observations with an orbital model 
when the astrometric wobble caused by the companion is below the noise level can have rather unexpected consequences. With 
new astrometric missions coming out within the next ten years, it is worth investigating the orbit reconstruction capabilities of 
such instruments at low signal-to-noise ratio. This is especially important because some of them will have no input catalogue 
thus meaning that all the orbital parameters will have to be derived from scratch. The puzzling case of almost parabolic orbits 
is also investigated. 
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1. Introduction 

Among the 23882 binaries of the Hipparcos Catalogue (ESA 
1997), only 235 were processed with an orbital model (with 
seven additional parameters), the so-called DMSA/O entries. 
However, the seven parameters were all fitted for only 45 of 
them. In the other cases, the value of some parameters were as- 
sumed from previous investigations, mainly spectroscopic and 
interferometric orbits. If the same proportion is assumed for 
GAIA, one will end up with about 5 10 5 systems for which 
one needs to derive an orbit from scratch. However, simulations 
have shown that the number of systems obser ved by GAI A for 
which an orbit is worth deriving is about 10 7 ( ESA 2000 ). 

The significant improvement of the fit of DMSA/O obser- 
vations when the orbital model is adopted does not mean they 
can constrain all the seven parameters. Assuming the value of 
some parameters was not just useful, it was sometime neces- 
sary in order to derive a realistic orbit. For instance, HIP 85749 
is a well known spectroscopic binary ( ^ucke & Mayor 1982 ) 
with a period of 418 days. Though almost three orbital rev- 
olutions took place during the Hipparcos mission and an or- 
bital model is well appropriate, a 7 -parameter \ 2 fit leads to 
e « 0.99 whereas the radial velocities yield 0.21. 

We report on the robustness of the fit of the seven parame- 
ters at different levels of noise on two dimensional and one di- 
mensional (Hipparcos-like) observations. From now on, it will 
be assumed that these observations have already been corrected 
for the parallactic effect and the proper motion, i.e. one will 
only deal with relative positions whose origin (center of mass 



or primary) is at rest. Unlike Sozzetti et al. (2001), we limit 
ourselves to S/N ~ 1. 

2. Synthetic material 

One thousand samples of No observations (No — 60 or 150) 
are generated by adding a Gaussian noise to the synthetic data 
of HIP 85749 (Table Q). The photocentric orbit is assumed to 
correspond to the astrometric one (ao = a\ = a). The (groups 
of) observations (Fig. [j]) are uniformly distributed over a pe- 
riod of Nu days (No = 1000 or 1600). Two scanning laws 
are also investigated. In one case, the orientation is uniform 
over 27T and the observations are uniformly distributed over No 
days. In the second, the system is observed on three consecu- 
tive revolutions of the instrument (the revolution period is set 
to 0.04 days), with the scanning direction slightly shifted (0.08 
radian). The groups of three observations are then uniformly 
distributed over Nd days. These properties, although realistic, 
do not specifically correspond to a particular forthcoming mis- 
sion (DIVA, FAME, GAIA & SIM). 

The transformation of the Cartesian coordinates (x, y) into 
the rectangular system attached to the instrument is: 



v = x cos t + y sin r 
w = — x sin r + y cos r 



(1) 

(2) 
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Several levels of noise (<j v — a w ) are investigated also: 0.25, 
0.5, 1, 1.5, 2.0 and 2.5 mas. Even if these values may sound 
very large with respect to modern instruments, they should be 
compared to the size of the projected orbit or a. From now on, 
2D and ID solutions will refer to orbits derived from (v, w) 
and v only respectively. 
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Fig. 1. The scanning direction makes an angle of r with respect 
to the x axis, v is the signed projection along the scanning di- 
rection of the separation between the two components (or the 
component and the center of mass), w is the projection along 
the cross-scan direction. 



Table 1. Orbit used to generate the observations. ABFG are the 



Thiele-Innes constants (Heintz 1978) 



Element 


Value 


A (mas) 


+8.666885e-01 


B (mas) 


-1.377152e+00 


F (mas) 


+3.176942e+00 


G (mas) 


-7.202854e-02 


e 


0.21 


P (days) 


418.42 


T(JD) 


2443723.5 


a (mas) 


3.3264 



For each pair of observations (ID and 2D), two least-square 
orbits are fitted using a local search only. The starting point is 
always the orbit given in Table |lj i.e. the orbit used to gen- 
erate the data. This is an unrealistic situation where the true 
orbit is known and we expect it to be the solution of the local 
search procedure. The distance between the solution of the lo- 
cal search and the true orbit is just a rough optimistic estimate 
of the reliability of the method. Indeed, a global search could 
find a better solution (in terms of \ 2 ) further away. However, 
if one does already run into troubles with the local search, one 
can seriously start worrying for the future. 

3. Orbit reconstruction 

3.1. 2D vs 1D data 

The distribution of P, e, and a are given in Fig. |2|. The dete- 
rioration of the distribution from 2D to ID and for decreasing 
S/N is quite obvious. For instance, the interval of periods with 
ID at 1 mas is about as large as with 2D at 1.5 mas. So, we do 
loose precision by only keeping the project ion along the s can- 
ning direction. However, as pointed out by Monet (1979 ), the 
period remains rather well defined since it is independent of the 
Keplerian nature of the signal. 




Fig. 2. Distribution of the derived e, P and a from 2D (thin 
line) and ID (thick line) data for different levels of noise with 
No = 1600d, No = 150 and the observations are grouped by 
3 

The behavior of the eccentricity is much worse. Even if ID 
and 2D distributions are still symmetrical at 1.5 mas, it is no 
longer the case at 2 mas and they start spreading all over the 
domain. However, instead of becoming uniform over [0, 1[, one 
notices an accumulation of circular and nearly parabolic orbits 
(lower left corner of Fig. ||). As already mentioned, there is no 
way a global optimization method can improve such a behavior. 

Even if e gets close to 1, an adequate choice of i and oj 
could nevertheless lead to a reasonable apparent orbit. For in- 
stance, a very eccentric orbit with a large semi-major axis will 
look likely if lj and i are close to | . The lower right panel of 
Fig. ^ exhibits such orbits with large semi-major axis. 

The results obtained from the different combinations of 
Np, Nq, scanning laws and noise levels are summarized in 
Table ||. As expected, the width of the distributions of e, P and 
a increases with the observation uncertainty. Although the re- 
lation between the width and the uncertainty is rather smooth, 
almost linear, there is a noise threshold above which the dis- 
tribution width, i.e. formally the parameter uncertainty, does 
explode. In the smooth regime, the improvement caused by the 
second coordinate ranges from 40 to 60%. The effect of the 
scanning law is barely noticeable. 

The period covered by the observations and their number 
do affect the width of the distribution more strongly. For in- 
stance, one remains in the smooth regime for all the 2D orbits 
with No = 150. As expected, the uncertainty on P decreases 
with the length of the observation campaign. However, the dis- 
persion of the observations seems to play a role also, especially 
in the ID case. Thus, 150 observations uniformly spread over 
1600 days give a better result, in terms of the region where the 
behavior of the distributions is smooth, than the same number 
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Table 2. Dispersion of the derived e, P, & a with respect to the number of observations (No), the duration of the observation 
campaign (N& days) and the scanning law (SL). SL=0 means that the No points are uniformly distributed over Npj. SL=1 
corresponds to observations grouped by three, the groups being uniformly distributed overs Nd- Regardless of SL, the scanning 
direction is uniformly distributed over 2n. 
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0.25 
0.50 
1.00 
1.50 
2.00 
2.50 
0.25 
0.50 
1.00 
1.50 
2.00 
2.50 
0.25 
0.50 
1.00 
1.50 
2.00 
2.50 
0.25 
0.50 
1.00 
1.50 
2.00 
2.50 



0.01 
0.03 
0.06 
0.09 
0.13 
0.20 
0.01 
0.02 
0.04 
0.06 
0.08 
0.10 
0.02 
0.03 
0.06 
0.09 
0.14 
0.18 
0.01 
0.02 
0.04 
0.06 
0.08 
0.10 



1.31 
2.65 
5.44 
7.84 
10.85 
14.29 
0.83 
1.62 
3.25 
5.15 
6.70 
8.29 
1.25 
2.58 
5.24 
7.62 
10.96 
13.54 
0.82 
1.66 
3.26 
4.95 
6.90 
8.36 



0.05 
0.10 
0.20 
0.33 
0.48 
777.02 
0.03 
0.07 
0.12 
0.19 
0.27 
0.33 
0.05 
0.10 
0.21 
0.31 
164.26 
615.85 
0.03 
0.06 
0.13 
0.19 
0.27 
0.34 



0.02 
0.05 
0.09 
0.15 
0.23 
0.31 
0.01 
0.03 
0.05 
0.09 
0.11 
0.17 
0.02 
0.05 
0.10 
0.17 
0.28 
0.32 
0.01 
0.03 
0.06 
0.09 
0.12 
0.18 



1.86 
3.91 
7.84 
12.01 
16.38 
25.39 
1.17 
2.37 
4.67 
7.11 
9.62 
12.02 
2.03 
4.21 
8.41 
13.22 
20.81 
26.13 
1.16 
2.40 
4.83 
7.35 
10.55 
13.19 



0.07 
0.15 
0.29 
364.70 
1266.08 
2350.59 
0.05 
0.09 
0.18 
0.29 
0.39 
650.23 
0.08 
0.17 
23.84 
536.81 
1443.31 
2062.24 
0.05 
0.09 
0.19 
0.28 
93.48 
743.73 



0.01 
0.03 
0.06 
0.09 
0.12 
0.18 
0.01 
0.02 
0.04 
0.06 
0.07 
0.10 
0.01 
0.03 
0.06 
0.09 
0.13 
0.17 
0.01 
0.02 
0.04 
0.05 
0.08 
0.09 



0.75 
1.50 
3.20 
4.68 
6.30 
8.37 
0.48 
0.95 
1.88 
3.01 
3.78 
4.79 
0.75 
1.50 
3.01 
4.77 
6.67 
8.26 
0.48 
0.92 
2.00 
2.99 
3.88 
4.91 



0.05 
0.10 
0.20 
0.30 
0.44 

624.03 
0.03 
0.06 
0.12 
0.18 
0.24 
0.33 
0.05 
0.10 
0.20 
0.30 
4.20 

429.38 
0.03 
0.06 
0.12 
0.19 
0.25 
0.31 



0.02 
0.04 
0.08 
0.14 
0.25 
0.30 
0.01 
0.03 
0.05 
0.08 
0.11 
0.16 
0.02 
0.05 
0.09 
0.18 
0.25 
0.29 
0.01 
0.03 
0.05 
0.07 
0.12 
0.17 



1.12 
2.22 
4.63 
7.32 
10.22 
13.29 
0.69 
1.36 
2.86 
4.21 
5.57 
7.48 
1.25 
2.48 
5.12 
8.86 
11.75 
15.52 
0.72 
1.36 
2.91 
4.42 
5.93 
8.00 



0.07 
0.15 
0.28 

344.91 
1335.89 
2101.33 
0.04 
0.08 
0.17 
0.25 
0.37 

809.25 
0.07 
0.15 
0.32 

361.40 
1125.18 
1766.75 
0.04 
0.09 
0.18 
0.27 

224.68 

732.08 



of observations grouped by 3. Actually, the latter case almost 
corresponds to 50 'normal' 2D points uniformly spread over 
the same period. 

3.2. 'Continuous' transition from 2D to 1D 

So far, we have considered the ideal case of 2D observations 
with the same precision on the two axes in the one hand and 
the poor case of ID data in the other. We now consider 2D ob- 
servations with different precisions along the two axes (along 
and cross scan). Indeed, a slightly worse cross scan precision 
might still prevent from running into the same troubles as with 
ID data at 2. 5 -mas noise level. 

As seen in Table |^, the sigma of the distribution of P and 
a can literally explode. Moreover, o is not necessary a good 
indicator since there is no guarantee that the distribution is 
Gaussian, nor even symmetric. In order to investigate that de- 
crease of the cross-scan precision, we plot (Fig. |3|) the semi- 
interquartile (a more robust indicator than a) of the distribu- 
tion versus the noise on w (cross scan direction) reckoned in 
<r v . Only three different a v are adopted since the point is just to 
see whether the overall behavior of SIQ(a w /a v ) is indepen- 
dent of the noise or not. 

The first thing to notice is that the shape of SIQ(a w /a v ) is 
a flat line whose height corresponds to the ID value. Another 



important point is the very steep growth at small a w /a v . The 
transition between 2D and ID is thus very rapid. The standard 
deviations of the distribution of e, P, and a exhibit the same 
rapid change at small abscissae. 

In Table |[ we have shown that, regardless of the period 
duration, the number of observations, . . . there is a threshold on 
a above which the behavior of a a explodes. For Table |] (and 
Fig. Q), the investigation is limited to a 1600-day mission, with 
150 observations grouped by three and a smaller increment for 
a between 1.5 and 2.0 mas is adopted in order to improve the 
identification of that threshold. 

The uncertainty on e and P looks like a linear function of 
the noise. The same is true for a up to about 1.9 mas (S/N ~ 
1.75). However, the ultimate goal of orbit reconstruction is to 
derive masses. So, down to which S/N can one go and still 
derive a useful result? The relation between a (actually at), the 
mass of the two components (Mi and Ma), P (in year) and the 
parallax w (reckoned in the same unit as a) is given by: 

a = M 2 (Mi + M 2 )- 2/3 mP 2/3 . (3) 

Therefore, if w and P are assumed to be error-free, M2 (Mi + 
M 2 ) _2/3 has the same relative uncertainty as a. Hence, al- 
though the fitted orbit still makes sense at S/N ~ 1.9, the 
theoreticians interested in masses should impose a larger S/N 
in order to obtain useful results. 
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Table 3. Relative precision of e, P and a for some S/N levels 
(i.e. a/a) 



a / a 



Fig. 3. Evolution of the semi-interquartile of the distribution of 
e, P and a as a function of the cross-scan precision (reckoned in 
along-scan unit). Filled squares, filled and open triangles cor- 
respond to a v = 2.5, 2.0 and 1.0 mas respectively. The long- 
dashed, short-dashed and continuous lines indicate the corre- 
sponding uncertainty with the ID solution. 

What does this lead to with GAIA? Like Lattanzi et al. 
(2000), we assume a noise level of 10/ias for V < 12 mag 
and a threshold at 1.9 for S/N. We can thus expect to derive 
an orbit whose a is as small as 19/ias. The relation between P 
and A/2 for objects with masses between 10~ 4 and 80 Jupiter 
masses, orbiting a solar-type star is plotted in Fig. || Our solar 
system is also represented. Observed from a Cen, about 1.3 pc 
away, the four Jovian planets are detected but the periods are 
too long with respect to the mission duration. In the other hand, 
the four inner most planets are still too light weight to cause a 
noticeable astrometric wobble. 

4. Global vs local optimization 

Up to now, it was assumed that the actual orbit was available 
as an initial guess of the local minimization. How good is this 
initial guess? The simulations in Sect. || show that the original 
orbit seldom minimizes the residuals (it is not the minimum of 
the least-square problem). However, such simulations do not 
tell us whether that original solution is always in the vicinity 
of the lowest minimum or not. Is the solution obtained with 
the local search from the original orbit the best one in terms of 
least-squares? How reliable is the orbit with the lowest mini- 
mum? 

The seven orbital parameters (^4, B, F, G, e, P, T) can be 
split into three groups: (a) those appearing in a linear way in 
the x 2 expression, namely A,B,F& G; (b) P because it can 
be derived reg ardless of the the Keplerian nature of the problem 
(Monet 1979); (c) e and T. This distinction is therefore propa- 
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Fig. 4. Evolution of the standard deviation of the distribution of 
e, P, and a at several noise levels 

gated to the way the \ 2 is minimized. The following scheme is 
thus adopted: 

1. Guess P usin g a period search technique (e.g. H ome & 
Baliunas 1986); 

2. Explore the (e, T)-space and derive A, B, F, G as the 
unique minimizer of x 2 when e, P, T are fixed. 

For each trial of e and T, another value of \ 2 i s evaluated. In 
this approach, the dimension of the global optimization prob- 
lem reduces to two. That transformation speeds up the mini- 
mization process and increases the chance of getting the best 
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0.1 
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0.001 1 1 1 — 1 1 1 

0.0001 0.001 0.01 0.1 1 10 100 

M3 (mass of Jupiter) 

Fig. 5. Limits of orbit reconstruction for objects orbiting a 
solar-type star at 250 pc (solid line), 100 pc (long-dashed line) 
and 1.3 pc (short-dashed line) with a = 19pas 

minimum. Several sophisticate methods exist for global search 
but, for the present purpose, a uniform 100 by 25 grid over 
[0, l[x [to, to + P[ where to is the time of the first observation 
seems to do just fine. The orbit thus derived can then be used 
as the initial guess of a local search in the seven dimensional 
space of the parameters. 

Once again, the method is assessed on synthetic observa- 
tions. However, contrary to what happened in Sect. |J we limit 
ourselves on realistic cases that are likely to be difficult: an ob- 
servation campaign of 1600 days, 150 observations grouped by 
3 and a noise of 2.5 mas. 1000 such sets of data are again gen- 
erated. For the sake of simplicity, we assume that the period 
search procedure returns the actual period (even if the 7p local 
search will update the period estimate later on). The minima 
obtained from the true orbit in the one hand and from the grid 
search in the other were then derived. The assumption about 
the efficiency of the period search technique might sound a bit 
optimistic and, indeed, it is. However, as we are going to see, 
even if such an unexpectedly good estimate of P is available, 
one can still run into troubles with other parameters. 

The two minima are different in 22 cases only. In two of 
them, the grid search gives a worse minimum than the one 
derived from the true orbit, thus indicating that a denser grid 
would give a better result. However, with the present sparse 
grid, one already notices that 2% of the solutions are not in the 
vicinity of the actual orbit. A denser grid can only increase that 
percentage. 

5. Almost parabolic orbits 

As already pointed out in Sect. |3.lj periodic though almost 
parabolic orbits (i.e. e w 1 and large a both compensated by 



i as 7r/2) show up quite often even when the assumed or fitted 
orbital period does not exceed the mission duration. When the 
seven orbital parameters are fitted, such orbits are quite usual 
and this is annoying and puzzling. For instance, it is annoying 
because that would strongly bias the distribution of e towards 
1 if such solutions were published. In the other hand, there 
seems to have no argument but statistical ones to discard these 
parabolic solutions. Indeed, the apparent orbit looks realistic 
(the high inclination and eccentricity cancel out). For instance, 
the two plots in Fig. ^| look similar although the eccentricity 
in the right panel is 0.998 where it is only 0.2 (spectroscopic 
value) in the left one. 

Fitting the Hipparcos Intermediate Astrometric Data of the 
genuine single star HIP 1 15331 with a 12-parameter model (i.e. 
the five astrometric parameters plus the orbital ones) leads to 
e = 0.9999. However, in that case, the Thiele-Innes constants 
are not significantly departing from (Pourbaix & Arenou 
2001). Thus, even if one almost derives a parabola, the whole 
solution is discarded so it does not matter. In the other hand, 
there are binaries for which the Hipparcos data leads to likely 
solutions (consistent with the spectroscopic orbits) even when 
the 12 parameters are adjusted from scratch (e.g., HIP 36377). 

With Hipparcos, besides the single stars and the binaries 
where the whole set of orbital parameters makes sense, there 
are many spectroscopic binaries for which an orbital model 
is worth adopting but the eccentricity must be assumed other- 
wise it gets close to a parabolic solution. Here, unlike the single 
star case, the Thiele-Innes constants are significantly non-zero. 
Actually, when the seven parameters are fitted, F and G get 
much larger than when e, P and T are assumed. In order to in- 
vestigate the reasons of such a strange behavior, let us remind 
the equations of motion in the projected orbit ( frleintz 1978 ): 

x = A(cos E - e) + Fy/ 1 - e 2 sin E, (4) 

y = B(cosE-e) +Gy/l - e 2 sin E (5) 
where E is solution of 

E-esmE = 2TT/P(t-T) (6) 

From these equations, it is clear that when e gets close to 1, F 
and G are no longer constrained. Thus, e ks 1 is a convenient 
way of getting rid of two parameters. Is e w la cause or a 
consequence of large F and Gl 

The relation between the projected areal constant V and the 
areal constant in the true orbit V is given by 

2r' = p 2 d0/dt = {AG - BF)y/l- e 2 2it/P 

= a 2 cosi\/l -e 2 2n/P (7) 
= r 2 dv/dt cos i = 2T cos i 

where v is the true anomaly in the true orbit. It is actually be- 
cause the law of areas holds in both the projected and true orbits 
that, in theory, the inclination can be derived. 

With precise 2D observations of rather short period bina- 
ries, r' is well constrained and so is cos i (the other elements 
being derived regardless of the law of areas). When S/N de- 
creases, r' gets less constrained. If the observations are clus- 
tered, r' is, on average, close to which leads to a parabolic 
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-35.6mas ^=262.54184331 deg. 35.6mas -42.3 mas ^*=262.54 1 84355 deg. 42.3ma S 

Fig. 6. An orbital model is worth adopting for HIP 85749 but rather similar apparent orbits nevertheless yield very different true 
orbits (the eccentricity is 0.2 and 0.998 for the left and right panels respectively). 



or edge-on orbit. As already mentioned, the former has the ad- 
vantage of getting rid of F and G in once. 

With ID data, the situation is worse. Indeed, there is no way 
to estimate T' and, therefore, i cannot be derived unless: 

- e is set (for instance, thanks to a spectroscopic estimate) 

- or successive ID data obtained with quite different scan- 
ning directions constrained 2D positions. 

The latter case happens when, for instance, the period is long 
enough for the actual position of the companion to remain un- 
changed in between distant observations and short enough for 
the orbit to be well covered during the mission. One way to 
achieve that would be to shorten the precession period of the 
satellite such that two consecutive observations would have 
quite distinct scanning directions. Unfortunately, the preces- 
sion period is no longer a free parameter in the design of the 
forthcoming missions. Another favorable circumstance is when 
pairs of distant observations are separated by a multiple of the 
period. This is rather unlikely but so seems to be the chance of 
fitting the Hipparcos data with a twelve parameter model from 
scratch. 

In the case of GAIA, its spectroscopic capability should 
help constraining the eccentricity of, at least, the short period 
binaries. Unfortunately, owing to the precision, this capability 
will be almost useless for the eccentricity of extra-solar planets. 

6. Conclusions 

One cannot expect to do as much and as well with ID obser- 
vations as with 2D ones. In between these two extremes, the 
second coordinate with a lower precision is not worth getting 
unless the ratio of the two precisions does not exceed 5. 



The x 2 minimization can advantageously be replaced with 
a 3-stage evaluation thus yielding a global search in a 2- 
dimensional (e, T) space. The period can indeed be indepen- 
dently guessed using period-search technique. Owing to this 
low dimension, a grid approach is shown to be quite efficient 
even if one cannot prevent some false solutions to show up. 



The results by Han et al. (2001 ) have already warned the 



community about the reliability of the astrometric orbits, es- 
pecially the inclinations, at low S/N even when some orbital 
parameters are adopted f rom spectroscopy ( pourbaix 200 1 ; 
pourbaix & Arenou 2001 ). The situation is likely to get worse 
when all the orbital parameters are derived from scratch. All 
these improvements in the way the orbits are derived will re- 
main useless unless one first finds a criterion for assessing the 
actual constraint on the projected areal constant and therefore 
on i and e. 
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